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ABSTRACT 

We search for z ~ 10 galaxies over ~160 arcmin^ of WFC3/IR data in the Chandra Deep Field 
South, using the public HUDF09, ERS, and CANDELS surveys, that reach to 5cr depths ranging from 
26.9 to 29.4 in i^ieo AB mag. z > 9.5 galaxy candidates are identified via J125 — i^ieo > 1.2 colors 
and non-detections in any band blueward of J125. Spitzer IRAC photometry is key for separating the 
genuine high-z candidates from intermediate redshift {z ~ 2 — 4) galaxies with evolved or heavily dust 
obscured stellar populations. After removing 16 sources of intermediate brightness (i^ieo ~ 24 — 26 
mag) with strong IRAC detections, we only find one plausible z ~ 10 galaxy candidate in the whole 
data set, previously reported in Bouwens et al. (2011). The newer data cover a 3x larger area 
and provide much stronger constraints on the evolution of the UV luminosity function (LF). If the 
evolution of the z ~ 4 — 8 LFs is extrapolated io z 10, six z ~ 10 galaxies are expected in our data. 
The detection of only one source suggests that the UV LF evolves at an accelerated rate before z ~ 8. 
The luminosity density is found to increase by more than an order of magnitude in only 170 Myr from 
z~10toz~8. This increase is > 4x larger than expected from the lower redshift extrapolation of 
the UV LF. We are thus likely witnessing the first rapid build-up of galaxies in the heart of cosmic 
reionization. Future deep HST WFC3/IR data, reaching to well beyond 29 mag, can enable a more 
robust quantification of the accelerated evolution around z ~ 10. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: luminosity function 



; 1. INTRODUCTION 

In recent years, great progress has been made in 
' quantifying the evolution of the galaxy population at 
the end of cosmic reionization around z ~ 6. Deep 
Hubble Space Telescope (HST) legacy fields, such as 

■ the Hubble Ultra -Deep Field (HUDF; j Beckmthelay 
! 12006) or GOODS (Gi avalisco et al.ll200 4l and wide area 

ground-based imaging, have made it possible to study 
the evolution of the UV luminosity functi on (LF) across 
z ^ 4 — 6 t o great accuracy (e.^. iBou wens et al.l 120071: 

■ lOesch et all |2QQ7|; IfeataeL^ UDM lO uchi et al. 2004; 
; iSawicki fc Th omDson"2006VMcLu re"eral.ii2009i ). 

I Over the last two years, with the installation of 
the Wide Field Camera 3 (WFC3) onboard the HST, 
the observational frontier of galaxies has now been 
pushed into the reionization epoch, as deep WFC3/IR 
data led to the identification of more t han 100 galaxy 
candidates at 2; ^ 6 . 5 — 8. 5 Te.g. IBouwens et al.l 
201 Ibt „B ouwens et al. 2010a'; 'Oesch et al. 2010d: 
McLure et al. 201 ; Bun ke r^t al. 2010; F i nkelstei n et ajj 
201(1 lYan et al.l I20T0I : iWilkins et aO [2OIOL [20111 
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iLorenzoni et al.ll2011l : [Trenti et al.ll201lD . This is essen- 
tial for estimating the contribution of galaxies to cosmic 
reionization. One of the most important conclusions 
from these studies is thus the realization that the 
UV luminosity density (LD) emitted by the galaxy 
population gradually falls towards higher redshift s. For 
example, the LD of the z ^ 3 galaxy population is about 
an order of magnitude larger than that of the 2; ~ 8 
population, about 1.5 Gyr earlier. 

How this evolves to even higher redshifts is still very 
unclear. A sizable galaxy population at z > 9 is ex- 
pected based on the first estimates of stellar popula- 
tion ages of z ~ 7 — 8 galaxies, indicating that these 
sources ver y lik ely started forming star s already at z ^ 
10-12 (e.R.Eab be et al. 2010b,a; Finkel stein et al .1120101 : 
iGonzalez et al.l [2 01 Of ). This is still somewhat uncer- 
tain due to possible nebular line emission contaminat- 
ing t he Spitzer photometry (e.g. iSchaerer fc de BarrosI 
l2010f ). Nonetheless, an early epoch of star- formation is 
also required by the mean redshift o f reionization as mea - 
sured by WMAP {zr = 10.6 ± 1.2; Komatsu et al.ll201lD . 
if galaxies are assumed to be the main drivers for this 
process. 

However, previous searches for 2: > 9 sources in the 
pre-WFC3 era only resulted in very small samples of rela- 
tively low reliability candidates, none of whic h have been 
confirmed (e g;, [go uwens et al. 2005; Stark et al.l 120071 : 
iHenrv et al.l 12008. 2009; Richard et al. .20081 ). This is 
mainly due to the extreme faint ness of the z ^ 10 galaxy 
population. Not only are these galaxies fainter due to 
their increased distance, but also they are expected at 
intrinsically lower luminosities. Additionally, the detec- 
tion of such high-redshift sources is further complicated 
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by the fact that they are invisible in optical data. Due 
to the highly neutral inter-galactic medium before the 
end of cosmic reionization, their UV photons are ab- 
sorbed shortward of the redshifted Lya line, which shifts 
to > Ifim dit z >?. Thus, these sources can only be seen 
in the NIR, where previous detectors were significantly 
lagging behind optical technology. 

With 40 X higher efficiency relative to NICMOS to de- 
tect high redshift galaxies in the NIR, WFC3/IR has the 
potential to change this and to push galaxy studies to 
beyond z ^ 9. Several deep and wide area WFC3/IR 
data sets have been taken already and several more are 
upcoming. The challenge for identifying genuine z > 9 
sources in these data sets, is that these galaxies will only 
be visible in one band (Hiqq). This has already led to 
some controversy in the ffist searches for z ^ 10 sources 
in the ffist-epoch WFC3/IR data ove r the Hubble Ul- 
tra D eep Field (see e.g. Bouwens et al.i l20l"Tal : lYan et aH 
l2010f ). In our recent analysis, which includes the full 
two-year WFC3/IR data over the HUDF as well as the 
shallower, but wider E arly Release Science (ERS) data, 
[Bouwens et alj (|2011af ) found only one single galaxy can- 
didate detected at > 5cr with an estimated redshift at 
z ~ 10.3. Given that about three should have been de- 
tected, if the evolution of the LF continued as extrapo- 
lated from the trends established across z ~ 4 to 2: ~ 6, 
this provided ffist tentative evidence for an accelerated 
evolution in the galaxy population from z ^ 8 to z ^ 10. 

In this paper we significantly expand on our first 
z ^ 10 analysis fro m WFC3/IR data presented in 
iBouwens et al.i ()2011a[ ) by extending the search to all the 
deep WFC3/IR fields in the Chandra Deep-Field South 
area that have since become available. The inclusion of 
the two deep HUDF09 parallel fields is especially useful, 
since both reach just '^0.5 mag shallower than the ultra- 
deep HUDF field but triple the search area for ^28 — 29 
AB mag sources. Additionally, we use different analysis 
tools developed by the first author that provide an inde- 
pendent analysis of the HUDF and ERS data. While the 
Lyma n-Break approach is similar in principle to that of 
IBouwens et al.l ((2011a) the use of independently tested 
software and procedures for the source detection and its 
analysis provides confirmation and validation. The ex- 
panded data set also covers > 3 x the area at moderately 
deep ~ 26.5 AB mag, thanks to the inclusion of the first 
epochs of CANDELS data over these fields. This will be 
used to constrain the evolution of the galaxy population 
over the ^ 200 Myr from z ^ 10 when the universe was 
- 500 Myr old to z - 8 at - 700 Myr. 

We start by describing the full data set in ^ and 
present the z ^ 10 candidate selection and its effi- 
ciency in ^ In ^ we present our new constraints 
on the LF at z - 10. We wih refer to the HST fil- 
ters F435W, F606W, F775W, F850LP, F098M, F105W, 
F125W, F160W as ^435, 1^606, ^775, ^850, ^098, ^105, -^125, 
i^i6o, respectively. Throughout this paper, we adopt 
Qm 0.3, Qa 0.7, Ho 70 kms~^Mpc~ \ i.e. h = 0.7. 
Magnitudes are given in the AB system (lOke fc GunnI 
[1983,). 

2. THE DATA 

Our analysis is based on the public WFC3/IR data 
sets that are availabl e over the GOODS South fields 
([Giavalisco et al.l 120041 ) as a result of three different pro- 




FiG. 1. The WFC3/IR fields (orange) around the GOODS 
South area used in this analysis (see 32 and Table [T}. The 
HUDF/HUDF09-1/HUDF09-2 fields comprise the deepest NIR 
data available to date, reaching to > 29 AB mag. These three 
WFC3/IR pointings are covered in a total of 192 HST orbits, which 
obtained photometry in Vios, J125, and ifieo- The ERS field con- 
sists of 60 orbits of WFC3/IR imaging in ^098, -^^125, and ifieo 
spread over a 2x5 tile. The CANDELS program is an on-going 
multi-cycle treasury survey consisting of two parts: CANDELS- 
Deep 60 arcmin^, hmited at Hiqq^ab ~ 27.5) and CANDELS- 
Wide (~ 40 arcmin^, limited at ifi60,AB ~ 26.9). Only J125 and 
if 160 data obtained before August 6, 2011 is used in this analysis 
(together with the shallow Yios imaging of CANDELS Wide). 

grams (HUDF09, ERS, and CANDELS), which we de- 
scribe below. The outhne of ah fields is shown in Fig- 
ure [TJ The key feature of these fields is that they have 
WFC3/IR coverage in J125 and i^ieo, which will be used 
to select z ^ 10 sources. Additionally, they have deep 
optical coverage as well as deep IRAC imaging, which 
is essential to exclude low redshift contamination (see 
section [3]). 

All the WFC3/IR data has been reduced follo wing 
standard procedures outlined, e.g., in Bouwens et al.l 
(|2011bl ). In particular, our reduction pipeline includes 
the subtraction of a super median image, careful image 
registration to the ACS frames and an automatic elimi- 
nation of pixels affected by persist ance. The final pixel 
scale of the images in our analysis is set to OVO6. A sum- 
mary of the HST data used in this analysis can be found 
in Table m The final resolution of the WFC3/IR data is 
- O^'ie (FWHM), and - 0'!09 in the optical ACS data. 

The ACS data on the GOODS fields we used are the 
v2 reductions that are publicly available from MASlfl 
(M. Giavalisco and the GOODS Team, in prepara- 
tion)^ The Spitzer data are the GOODS IRAC images 
([Dickinson et al.ll2003l ) a s made publicly ava ilable by the 
SIMPLE team (see e.g. iDamen et al.ll2009f ). Addition- 
ally, we include the newly acquired Spitzer IRAC [3.6] 

^ http:/ /archive. stsci.edu/prepds/goods/ 
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TABLE 1 

Summary of 5cr Depths^ of Observational Data Used in our Analysis 



Field Area [arcmin^] B435 1^606 ^775 ^^sso ^098/i05 -^125 ^160 



HUDF09 


4.7 


29.2 


29.6 


29.4 


28.8 


29.1 


29.3 


29.4 


HUDF09-1^ 


4.7 




28.9 


28.7 


28.6 


28.6 


28.8 


28.6 


HUDF09-2b 


4.7 


28.8 


29.3 


28.9 


28.7 


28.6 


28.9 


28.9 


ERS 


41.3 


27.8 


28.0 


27.5 


27.2 


27.4 


27.8 


27.6 


CANDELS-Deep^ 


63.1 


27.8 


28.0 


27.5 


27.2 




27.7 


27.5 


CANDELS-Wide^ 


41.9 


27.8 


28.0 


27.5 


27.2 


27.1 


27.2 


26.9 



^ Measured in circular apertures of 0^25 radius. 

^ New relative to IBouwens et al.l (|2011al ) for 2; ~ 10 search. 



and [4.5] data over the HUDF field from the lUDFlO sur- 
vey (proposal 70145, PI: Labbe). This so far adds ~ 130h 
of observations, which increases the depth in both filters 
by an additional ~0.4 mag. Where available, we also 
matched our WFC3/IR sources with the publicly avail- 
able GOODS - MUS IC multi-band photometry catalog of 
ISantini et aP (j2QQ9[ ). 

2.1. Full HUDF 09 Data Set 

The HUDFQ9 program (PI: Illingworth:' Bouwens et ahl 
l2Qllb[ ) consists of 192 HST orbits to provide ultra-deep 
WFC3/IR imaging over three pointings centered on the 
HUDF (Beckwith et aIll2QQ6[ ) and its two p arallel fields 
from the UDF05 program (PI: Stiavelli; lOesch et al.l 
120071 ). The program has been completed, providing the 
deepest IR images ever taken. It comprises 3 x 4.7 
arcmin^ imaging in the three filters Yios, ^125, and i^ieo, 
reaching down to Hiqq^ab = 29.4, 28.6, and 28.9 (5cr in 
0'.'5 diameter apertures) for the HUDF, HUDF09-1, and 
HUDF09-2, respectively (see Fig. [1]). For a more de- 
tailed de scription of this d ata set and the data reduction 
see IBouwens et al.l (j2011br ). 

2.2. Wide Area Data 

In addition to the ultra-deep HUDF09 data, we also 
analyzed shallower, wider area WFC3/IR imaging from 
the ERS and CANDELS programs, in order to constrain 
the volume density of more luminous star-forming galax- 
ies at z ~ 10. 

The Early Release Science Data (ERS) provide 
WFC3/IR imaging of 41 arcmin^ of the northern part 
of the GOODS South field. Two orbits of WFC3/IR 
imaging were obtained in each of the filters loos, ^125, 
and i?i60: over a 2 x 5 grid of pointings (60 orbits in 
total). These data are reduced in an analogous way 
to our HUDF09 data, and are aligned and drizzled to 
the GOODS ACS mosaics after rebinning to a OVO6 pix- 
elscale. These data reach to i^ieo = 27.6 (see also 
IBouwens et al.ll2011bi). For a more detaile d description 
of this data set see I Windhorst "eFall (|20TT1 ). 

The last two fields included in our analysis are ob- 
tained as part of the Mult i- Cycle Tre asury program 
CANDELS (PI: Fabe r /Ferguson; Grogi n et al.l [20TTI : 
iKoeke moer e t all 1201 IV In particular, we include the 
first six visits of the CANDELS-Deep program (obtained 
until August 6, 2011), which covers the central part of 
GOODS South in 3 X 5 tiles with ~ 6000 s exposures in 
both J125 and i^ieo in a total of 92 orbits. This data 
covers ~ 63 arcmin^ and reaches to i^i60,AB ~ 27.5 
mag. Additionally, we also included the imaging data 



of the supernova follow-up program of CANDELS (PI: 
Riess), which adds imaging over two pointings over CAN- 
DELS Deep (one of which is essentially centered on 
the HUDF). Finally, we made use of the 29 orbits of 
WFC3/IR data of the CANDELS- Wide survey (Yios, 
J125 and i^i6o, obtained until March 29, 2011). These 
comprise 9 WFC3/IR pointings (~ 42 arcmin^), complet- 
ing the coverage of the GOODS South field, and reach to 
Him,AB = 26.9 (- 2000 s exposures). As for the ERS, 
the WFC3/IR data of the CANDELS program has been 
aligned to the GOODS ACS mosaics with a pixel scale of 
0V06. The part of the CANDELS field overlapping with 
the WFC3/IR HUDF has been omitted when analyzing 
this data set in order not to duplicate the analysis of that 
area. 

We will subsequently refer to the combination of the 
ERS and the two CANDELS fields as 'Wide Fields'. 

3. SOURCE SELECTION 
3.1. Catalog Construction 

Sou rce catalogs are derived with the SExtractor pro- 
gram (jBertin fc ArnoutsI 119961 ). which is used to detect 
galaxies in the i^ieo images and perform matched aper- 
ture photometry on PSF- matched images. The colors 
used here are based on isophotal apertures derived from 
the HiQo images, and total magnitudes are measured in 
standard 2.5 Kron apertures, corrected by 0.2 mag in 
order to account for flux loss in the PSF wings. 

The detection significance of sources was established 
in OV25 radius apertures. The RMS maps were scaled 
based on the detected variance in 1000 random apertures 
for each WFC3/IR frame on empty sky regions after 3a 
clipping. This procedure ensures that the SExtractor 
weight maps correctly reproduce the actual noise in the 
images. Subsequently only sources with signal-to-noise 
ratios larger than 5 (in 0^25 radius apertures) in i^ieo 
are considered. 

3.2. Ji2b-dropout Candidate Selection 

Galaxies at z > 9.5 are expected to exhibit very red 
J125 — HiQo colors since the redshifted, strong Lya ab- 
sorption (by the predominantly neutral inter-galactic hy- 
drogen) cuts into the flux in the J125 filter (see Figure 
[2]). This makes such high redshift galaxies completely 
invisible blueward of J125, and we use this fact for their 
identification. 

As can be seen from Figure O however, also passively 
evolving or dusty galaxies at intermediate redshifts {z ~ 
2.5 — 4) can exhibit similarly red colors in J125 — i^ieo- 
While the requirement of optical non-detections removes 
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Fig. 2. — Color selection of z > 9.5 galaxies f ^3.2|) . Shown are 
the Ji25 — HiQo colors of different types of galaxies as a function of 
their redshift. Star-forming galaxies are shown as solid blue lines. 
The lighter blue spectral energy distribution (SED) is reddened 
by E(B-V)= 0.15 mag using a Calzetti et al. (2000) dust law. As 
the Lya absorption due to inter-galactic neutral hydrogen shifts 
into the J125 band, galaxies start to exhibit progressively redder 
colors beyond z ^ 8. A selection with J125 — i^ieo > 1-2 thus 
identifies galaxies at z > 9.5. Strong Balmer and 4000 A breaks 
in evolved 2; ~ 3 galaxies, combined with some dust obscuration 
can also result in very red J125 — Hiqq colors. The dashed lines 
correspond to more evolved galaxies dominat ed by progressively 
older stellar populations from the library of Coleman et al. ( 1980) 
as well as a 1 Gyr old single stellar population (SSP) from the 
library of Bruzual & Chariot ( 2003). Spitzer IRAC data provides 
a way to separate these different populations at intermediate and 
very high redshifts (see text and Fig. [S])- 
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Fig. 3. — Sensitivity limits (2a) of some of the data used in this 
analysis. The red arrows indicate the limits of the ACS, WFC3/IR 
and IRAC data over HUDF, while blue arrows correspond to the 
ERS data set. The SED of a 2; = 10 star-forming galaxy is plot- 
ted in purple, showing the distinguishing feature of complete ab- 
sorption shortward of the redshifted Lyo; line. Evolved and dusty 
galaxy SEDs can mimic the same features as a 2; ~ 10 source, 
i.e. red J125 — Hiqq and undetected in the optical data, unless 
this is extremely deep. Such sources are typically detected in the 
IRAC data. Alternatively, the solid and dashed green lines show 
possible SEDs of 2; ~ 2.7 galaxies which could escape the IRAC 
detection at 2cr. These correspond to evolved galaxies (500 Myr), 
with moderate amounts of dust extinction, and are normalized to 
the 5 cr detection limits of the ifieo bands. Note that much larger 
amounts of dust can be hidden for interloper galaxies at the de- 
tection limit of the HUDF (green solid. Ay = 2.0 mag) than for 
the ERS (green dashed. Ay = 0.6 mag). Thus even with IRAC it 
can be a challenge to separate real 2: ~ 9 — 10 sources from lower 
redshift contaminants. 



the bulk of lower redshift contamination, certain inter- 
mediate redshift galaxies with evolved or dusty stellar 
populations can still be included due to the fact that 
the optical data does not reach deep enough, if at sim- 
ilar depth as the IR (see Fig. [3]). Deep Spitzer IRAC 
data provides a way to identify contaminating galaxies. 
These are expected to exhibit very red i^ieo — [3.6/im] 
colors, which discriminates them from genuine z ^ 10 
candidates. 

To exclude possible low-redshift contamination, we 
thus use two steps to select z > 9.5 galaxy candidates. 
For the first step, the primary criteria are based on HST 
data only: 

(-^125 — ^16o) > 1-2 



S/N{Hieo)>5 A 5/7V(< J125) < 2 A 



X.opt ^ X.cut 

Additional to excluding objects that are detected in 
any band blueward of J125 at more than 2cr, we in- 
clude a cut in the optical Xoyt v^^l^^ of ^ galaxy (see 
e.g. iBouwens et al.llMTbl : iBouwens et al.ll2Qllal ). This 
is computed from the 0^.^25 radius aperture fluxes as 
Xopt — ^i^^^^ifi) {fi/^i)^i where the sum runs over all 
the bands available in the given data set blueward of 
J125, i.e. it includes all the available optical data as well 
as the NIR band Y105 for the HUDF09 data, and Yogs in 
the ERS. The relatively large apertures were chosen in 
order sample > 70% of the light of point-like sources. 

The limiting Xcn^ derived from photometric scat- 
ter simulations. They are set to exclude the majority 
of interlopers which remain undetected at 2cr purely due 
to photometric noise, but not to cut a substantial frac- 
tion of galaxies with real zero flux in the optical bands. 
The scatter simulations utilize all galaxies in our cata- 
logs that are 1 — 3 mag above the completeness limit, 
applying photometric Gaussian noise from 1 mag fainter 
sources. From these simulations it is clear that contam- 
ination is mainly an issue at 0.75 mag above the com- 
pleteness limits, but that 60 — 80% of contaminants 
can be eliminated by using a Xop^ hmit of Xcnt ~ 2.8 or 
2.4, for 5 fllters or 4 fllters, respectively. In the IIUDF09 
data, the resulting number of expected contaminants due 
to photometric scatter is thus reduced from ~ 0.5 source 
per WFC3/IR fleld to - 0.1 source. 

On the other hand, the adopted Xcnt liiiiits do remove 
an additional ~ 20% of sources with real zero flux, simply 
due to Gaussian statistics. This reduction of the real 
galaxy sample is reflected in our subsequent analysis in 
the reduction of the selection volume. 

All galaxies passing the above selection criteria, us- 
ing both the ACS and WFC3/IR data, are retained and 
analyzed individually. These total to 17 sources with 
Hi6o,AB in the range 23.6 — 28.8 mag; one source in the 
HUDF, none in the parallel HUDF09 flelds, three in the 
ERS and 8 and 5 in the CANDELS Deep and Wide, 
respectively (see Tables [2] and |4|) . 

Interestingly, only one source (the previously reported 
galaxy with ^leo ^ 29 mag from Bouwens et al. 2011) 
did pass our selection in the three deep HUDF09 flelds, 
while the shallower CANDELS and ERS flelds contribute 
a total of 16 sources (all with Hiqq^ab ^ 26 mag). 
Upon inspection of their images, it turns out that all 
these brighter sources are very well detected in the IRAC 
data, even in the shallow 8.0/im band. Their measured 
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TABLE 2 
The 2; ~ 10 Galaxy Candidate 



ID a 5 ifieo -^125 - ^i60 S/N/j^ieo 

HUDFj-39546284 03:32:39.54 -27:46:28.4 28.8 ± 0.2 > 2.02 6.3 



F435W F606W F775W F850LP F105W F125W F160W [3.6] [4.5] [5.8] [8.0] 
HUDFj-3954628r4 V ** K K % JpJt- >" ■ 

....... t 1 I V U-JsL^_ 'fl 

Fig. 4. — 5^^x5^^images of our only viable z ~ 10 galaxy candidate (see also lBouwens et al.ll20 11gJ). The images show, from left to right, 
^4355 ^606 5 ^7755 ^850, ^105, -^125, -^^160, IRAC [3.6], [4.5], [5.8], and [8.0]. The J125 and i/ieo images shown here are combined from both 
the HUDF09 and the CANDELS survey data (which add 0.05 mag in depth). The IRAC [3.6] and [4.5] images include the additional 
~ 130h of data obtained over the HUDF as part of the lUDFlO program so far. As can be seen, the source is only significantly detected 
in HiQQ. The IRAC [3.6/im] data shows a non-negligible excess of fiux in the vicinity of the candidate. However, this fiux is very likely to 
be associated with the brighter neighboring source to the upper right. After subtraction of all nearby sources in IRAC, the candidate is 
undetected in both IRAC bands at less than Icr. The images are oriented N top, E left. 



^160 — [3.6] colors are in the range 1.6 — 4.3, which, for 
a z ~ 10 source, would correspond to a UV continuum 
slope P > —0.2, or a dust reddening with Ay > 1.6 
mag. Given that galaxies at the bright end of the 
z ^ 7 population are measured to have very low ex- 
tinction values and continuum slopes of P — 2.0 ± 0.2 
(see e.g. B ouwens et al. 2010b; Finkelstein^eLalJ 120101 : 
iDunlop et al.ll201ll : IWilkins et al.ll2011 b). the extremely 
red colors of these galaxies rules out z ^9 solutions with 
any sensible SED. 

All sources with IRAC detections are thus removed 
from our sample of potential z ^ 10 galaxies, reducing 
the sample to one single candidat e in the HUDF, pre- 
viously reported in iBouwens et al.l ()2011ai see Figure H] 
and Table [2]). The 16 removed sources are shown in the 
appendix in Figure [TT] and listed in Table HI 

3.3. The z ^ 10 Galaxy Candidate 

Images of the only possible 2: ~ 10 galaxy candidate 
are shown in Figure HI The source is detected at 6.3 
a in HiQo (measured in circular apertures of 0'.'25 ra- 
dius). This is higher but completely consistent with the 
IBouwens et al.l (|2pila) significance estimates, which are 
based on smaller apertures. As can be seen, the source 
is not significantly detected in any other band. Its value 
of xlpt — 2.77 is very close, but just below the limit of 
xlut = 2.8. This is mainly due to a 1.5a fiux excess in 
^775, which appears to be due to an extended structure in 
the background of that image. When adopting smaller 
apertures, the Xopt value is found to be reduced, indicat- 
ing also that this excess of fiux is not associated with the 
source itself. 

The Spitzer IRAC 3.6/im data shows some fiux from 
a nearby source. However, after subtraction of all the 
neighboring sources in IRAC, the candidate is undetected 
(0.09<j), with a 2(j upper limit on its IRAC [3.6] magni- 
tude of > 27.2 mag AB (Gonzalez et al., in prep.). This 
thus corresponds to i^ieo — [3.6] < 1.6 at 2cr, which is 
much bluer (by > 0.4 mag) than the typical low-redshift 
contaminants that we culled from our sample (see also 
next section). After adding the newly acquired IRAC 
data from the lUDFlO program (Spitzer proposal 70145, 
PI: Labbe) to the GOODS IRAC data and removing 
neighboring sources, the z ^ 10 candidate is also un- 



detected at [4.5] {S/N = 0.3) providing added weight to 
the likelihood of it being at high rather than low reshift. 

We derive the photometric redshift of the candi- 
date using the co de ZEBRA (jFeldmann et al.l 120061 : 
lOesch et al.l 1201 Obf) with synthe ti c ste llar population 
models from iBruzual fc CharlotI (|2QQ1 ) to which we 
adde d nebular continuuni and line emission following, 
e.g. JSchaerer fc de BarrosI (|2009l ). Using the full 11 band 
fiuxes and fiux errors, we derive a photometric redshift 
for this source of Zphot — 

10.41q;4, with a likelihood of a 
low redshift solution at z < 8 of < 6%. The full spectral 
energy distribution (SED) of the source and its redshift 
likelihood function are shown in Figure [5l The best-fit 
SED corresponds to a very young, dust-free star-burst 
(see also Gonzalez et al. in prep.). 

The best low redshift solution is found at ziowz = 2.7, 
for an evolved, very low mass galaxy SED (M = 3 x 10^ 
Mq) with moderate extinction. Interestingly, this SED 
is expected to be detected only at ~ 1.5cr in J125, but it 
nevertheless has a significantly higher value (13.7 com- 
pared to Xbest ~ '^•^)- Deeper HST data shortward of the 
break would he extremely useful in order to constrain the 
possible non-detection of the source shortward of 1.4 /im. 
As with other high-redshift catalogs, the added shorter- 
wavelength optical/near-IR data would play a key role in 
helping to further tighten its photometric redshift mea- 
surement. 

3.4. The Dusty and Evolved Contaminants of 
Ji25-Dropout Selections 

The properties of the 16 intermediate brightness 
sources that did not pass the IRAC non-detection criteria 
are discussed in the appendix. From fitting their SEDs, 
these sources are found to be mostly massive galaxies 
(M > 5 X 10^° Mq) with obscured but evolved stellar 
populations at z 2 — 4. 

Interestingly, all these galaxies are essentially limited 
to HiQo ~ 24 — 26 mag. This is ~1 mag brighter than the 
detection limits of our bright survey fields (see Figure [10] 
in the appendix), and thus suggests that such very red 
galaxies have a somewhat peaked luminosity function. 
This will have to be confirmed and quantified with future 
wide area data. However, at face value, it would appear 
that contamination from such sources with J125 — ^leo > 
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Fig. 5. — Best-fit SED of the only viable 2; ~ 10 candidate, 
previously reported in Bouwens et al. ( 2011a). The measured pho- 
tometry is shown with the red circle with errorbars and with 2cr 
upper limits in case of non-detections. All fiuxes and fiux errors 
were used in the SED fit, however, even if they were negative. 
The best fit SED is found at Zp^ot = 10-4 and is shown as a solid 
line (Xhest ~ '^•^)' corresponding to a dust-free, young galaxy. The 
best template for a low redshift solution is also shown as the dashed 
green line. This SED has a redshift of ziowz = 2.7, and is passive 
with an age of 650 Myr and of rather low mass (only 3 x 10^ Mq). 
Additionally, the low-z SED is reddened with E(B — V) = 0.1 mag. 
Based on its larger value (x^owz ~ 1^.7), this SED is formally 
excluded at > 90% probability, however. The redshift probability 
function is shown in the inset in the upper left. The sharp decrease 
above z ~ 10.5 of p{z) is due to the use of a LP prior at 2; > 8 and 
due to the fact that Lya absorption starts affecting the ifieo band 
(which would require the source to be brighter intrinsically). 

1.2, and with i^ieo — [3.6] ^ 2 is less problematic at fainter 
magnitudes. We note, however, that contamination from 
similar galaxies with less extreme colors (which may be 
more abundant also at i^ieo > 26 mag) may still be non- 
negligible due to photometric scatter (see also ^33j). 

Furthermore, it is interesting to note that it is very 
difficult to construct clean z ~ 10 galaxy samples purely 
based on HST data alone. Even if we were to increase the 
color criteria to J125 — i^ieo > 2.0, there would still be 
two contaminating galaxies in the sample together with 
our only viable z ~ 10 candidate. Thus, also for future 
constraints on the bright end of the z ~ 10 LF, it will 
be important to perform additional follow-up studies to 
validate the candidates, e.g. with Spitzer. This will be 
less of a concern for future z ~ 9 galaxy samples, which 
can be obtained, e.g., based on new F140W filter data. 
In such a data set intrinsically red galaxies can be sorted 
out by requiring a blue continuum across F140W and 

^160- 

3.5. Possible Sources of Sample Contamination 

Here, we only give a brief summary of the possible 
sample contami nation. For a thorough discussion we re- 
fer the reader to iBouwens et al.l (|2011a[ ) . 

Essentially, the only probable chance for contamina- 
tion is due to photometric scatter of a red lower redshift 
source. We estimate this to be a 10% chance based on 
our photometric scatter experiments described in section 
13. 2| including the x^opt cuts. 

Other typical contaminants to LBG selections such as 
very cool dwarf stars and supernovae can essentially be 
excluded based on the relatively blue J125 — i^ieo (< 
1.1) colors of stellar SEDs, on the fact that the source 



is detected in both the first and the second year of the 
HUDF09 WFC3/IR data, and due to the fact that the 
candidate shows clear signs of an extended morphology. 

Additionally, it is very unlikely that this source is spu- 
rious, since the flux distribution in circular apertures ran- 
domly distributed over empty regions of the HUDF i^ieo 
image are nearly exactly Gaussian, and the source is well 
detected at 6.3cr. 

4. THE ABUNDANCE OF Z - 10 GALAXIES 

In this section we compute the expected abundance of 
z ~ 10 galaxies in our dataset, and derive constraints on 
the 2: ~ 10 LF based on our data. 

4.1. Completeness and Selection Functions 

In order to estimate the number of sources we expect 
in our data from a given LF, we have to estimate the 
completeness as a function of magnitude C{m) and se- 
lection function as a function o f redshift an d magnitude 
S{z, m). Following Oesch et al.l ([20071 12 0091 ). this is done 
by inserting artificial galaxies in the observational data 
and rerunning the source detection with the exact same 
setup as for the original catalogs. This is done for each 
of our fields individually. 

Two sets of simulat i ons w ere run. In the first set, we 
follow Bouwens et al.l ()2003f ) , where the artificial galax- 
ies are 'cloned' from the z ^ A dropout sample of the 
GOODS and HUDF fields. The images of these sources 
are adjusted for surface brightness dimming, the differ- 
ence in angular diameter distance, as well as a size scaling 
of (l + z)~^ as observed for the Lyn ian Break galaxy pop- 
ulation across z ^3 — 7 (see e.g. iFerguson eran[200l 
IBouwens et al.|[2004l : lOesch et al.ll2010a[ ). These are then 
inserted in the observed images with galaxy colors as 
expected for star-forming galaxies between z = 8 and 
z = 12. When computing the galaxy colors we assume a 
distributio n of UV continuum s lopes w i th /3 = — 2.5=b0.4 
/see e.g. ..Bouwens et all I2009L l2010bl : iFinkelstein et al.l 
I20T0I : iStanwav et al.ll2005[ [ 

From these simulations, we compute the completeness 
as a function of observed i^i60,AB magnitude for each 
field, taking into account the scatter and offsets between 
input and output magnitudes. Additionally, we com- 
pute the selection probabilities as a function of redshift 
and magnitude by measuring the fraction of sources that 
meet our selection criteria. By construction, galaxies are 
selected at > 50% at redshifts z > 9.5. 

In the second set of simulations, we repeat the above 
procedure. However, instead of using observed galaxy 
images, we use theoretical galaxy profiles from a mix of 
exponential (Sersic n = 1) and de Vaucouleur (Sersic 
n = 4) profiles. The size distribution is chosen to be log- 
normal, again with the same size scaling as a function 
of redshift. The completeness and selection functions 
of our two procedures are in excellent agreement. This 
demonstrates the reliability of our approach, which ap- 
pears to be essentially independent of the adopted galaxy 
profiles (unlike wh at has been claimed elsewhere, e.g. 
iGrazian et al.|[2QTQl but see also Bouwens et al. 2010b). 

4.2. The Expected Number of z ^ 10 Sources 

The expected number of sources in a given magni- 
tude bin nii can be estimated for any given LF, (l){M), 
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Fig. 6. — Expected number of 2; ~ 10 candidates in the dif- 
ferent fields assuming the LF to be as measured at 2; ~ 8 (top; 
[Bouwens et al.h2011b>) or as expected from extrapolating the z ~ 
4 — 8 trends to 2; ~ 10 (bottom). The one detected source is in 
stark contrast to the ~ 27 sources expected if the LF was constant 
across the 170 Myr from 2; ~ 10 to 2; ~ 8. The LF appears to 
evolve at a significantly accelerated pace with respect to the em- 
pirical evolution observed across 2; ~ 4 — 8. 6 sources at 2; ~ 10 
should be detected when extrapolating the lower redshift trends, 
while only one probable candidate is found. 

through: 

Ar;^P= / dm [ dz^S{m,z)C{m)(l){M[m,z]) 

J Am J 

The result of this calculation is shown in Figure [6] for 
ea ch individual field usi ng the observed LFs at z ~ 
8 (|Bouwens et al.l l2Qllbl ) , as well as that expected at 
z ~ 10. The latter is based on extrapolating the evo- 
lution of the Schechte r function parameters presented in 
iBouwens et al.l (|2Qllbi) , which is based on fitting the like- 
lihood contours of LBG LFs across z ~ 4 — 8. In particu- 
lar, significant evolution is only seen in the characteristic 
luminosity M*, which is found to dim by 0.33 mag per 
unit redshift. The other two Schechter function param- 
eters were kept constant, in agreement with the small 
ev olutionary trends that were found at low significance 
in IBouwens et al.l (|2011b[ ). For reference, the final pa- 
rameter evolution we use in this work is: 

M^{z) = -21.02 mag + - 3.8) x 0.33 mag 

(f)^ = 1.14 X 10~^ Mpc~^mag~-^ = const 

a = — 1.73 = const 

As can be seen in Figure [6l the wide- area, shallow data 
of CANDELS and ERS are very useful for constraining 
the evolution in the bright end of the LF. In particular, 
if there was no evolution in the LF from z~10toz~8. 



the wide area data should contain 8.5 z ^ 10 sources. In 
total, we would expect to see ~ 27 sources if the LF was 
unchanged over the 170 Myr from z ^ S to z ^ 10. This 
implies that the wide area data provides essentially 30% 
of the total search power in case of an LF evolution in 
^* only. 

Given that we only detect one probable candidate, the 
LF appears to drop faster than expected from the em- 
pirical lower redshift extrapolation. In particular, we do 
not detect the ~ 6 galaxies that we would have expected 
to find at z ~ 10 if the lower redshift trends remained 
valid. From these extrapolations we predicted to find 
three sources in the HUDF (cons istent with the expecta- 
tions from Bouwen s et al.ll2Qlla[ ) and about one in each 
of the two HUDF09 parallel fields. 

The Poissonian probability to find < 1 source, given 
that 6 are expected is < 2%. This remains significant, 
even after including cosmic variance, which adds an ad- 
ditional uncertainty on these low number counts of about 
45 — 50% for an individ ual WFC3/IR po inting (see e.g., 
iTrenti Stiave lli 2008; iRobertsod [2010l ) . We derive an 
upper limit of 6% to the probability of finding < 1 source 
in our search area, based on t he cosmic variance calcula- 
tor of iTrenti fc StiavelH ()2008l ) and combining the num- 
ber counts uncertainty in the different fields assuming 
the final distribution is Gaussian (justified by the central 
limit theorem). Therefore, the detection of an acceler- 
ated evolution relative to the low-redshift extrapolation 
is significant at > 94%. 

4.3. Constraints on the z ^ 10 Luminosity Function 

The accelerated evolution in the UV LF can also be 
seen from our constraints on the z ^ 10 LF. The step- 
wise LF is computed using an approximation of the ef- 
fective selection volume as a function of observed mag- 
nitude Feff(m) = dz^S{z,m)C{m). The LF in 
bins of absolute magnitude is then given by (j){Mi)dM = 
N^^^/Vefi{mi). This is shown in Figure [71 where the LF 
was evaluated in bins of 0.5 mag, and non-detections cor- 
respond to la upper limits including the effects of 50% 
cosmic variance per pointing. 

As can be seen in Figure [71 the current constraints 
on the 2: ~ 10 LF are significant even at bright mag- 
nitudes where the wide area data are particularly valu- 
able. These data sets reduce the upper limits by more 
than an order of magnitude relative to using only the 
HUDF09 fields, therefore indicating that the z - 10 LF 
at Mi4oo < — 20 drops by a factor ~ 4 — 5 with respect to 
the observed LF at z ~ 8. However, from these shallow 
data sets, no constraints can be obtained on any acceler- 
ated evolution of the UV LF around z ^ 8 — 10. This only 
becomes apparent at M1400 > —20 (corresponding to 
H16O, AB ^ 27.5 mag). Such faint limits are only probed 
by the HUDF09 data set. In particular, at M1400 = — 19, 
the upper limit on the LF is a factor ^ 3 below the ex- 
pectation (Fig. [7|). Thus, it is clear that data reaching 
to deeper than Hiqq^ab = 28.5 mag will be necessary to 
further constrain the drop in the LF in the future (which 
is beyond the reach of the current MCT programs). 

In order to quantify the change in the LF from z ^ 10 
to z ^ 8 more robustly, we consider two possible scenar- 
ios. First, we assume the accelerated LF evolution occurs 
only in M*, at a constant rate since z ~ 6, and we fit the 



8 



Oesch et al. 




Fig. 7. — Constraints on the z ~ 10 LF from our combined data 
set, evaluated in 0.5 mag bins. The upper hmits correspond to la 
Poissonian hmits including the additional uncertainty of 50% cos- 
mic variance per pointing. It is clear that the luminosity function 
evolves strongly from z ~ 8 to 2; ~ 10, as our upper limits are a 
factor ~2-5 below the measured 2; ~ 8 LF of Bouwens et al. (or- 
ange data; 2011b). The expected 2; ~ 10 LF as extrapolated from 
fits to lower redshift LBG LFs is shown as a dashed red line. Using 
this LF, we would expect to detect six sources in the full data set 
3.4 in the HUDF, and ~ 1.2 in each of the HUDF09 paral- 
lels; see Figure O. For comparison also the z ^ 4 and z ~ 6 LFs 
are plotted (Bou wens et al. 2007), showing the dramatic build-up 
of UV luminosity across ^ 1 Gyr of cosmic time. The light gray 
vectors along the lower axis indicate the range of luminosities over 
which the different data sets dominate the 2; ~ 10 LF constraints. 

Poissonian likelihood for the observed number of sources. 
This can be written as £ = JJ. Y[. P(A/'°^', A^,^^), where 
j runs over all fields, and i runs over the different mag- 
nitude bins, and P is the poissonian probability. 

Our extrapolation o f the UV LF is a mod ification of 
the fitting formulae of iBouwens et al.l (|2Qllbl ). We thus 
use = 1.14 X 10"^ Mpc~^ =const. and a{z) = 

— 1.73 =const., and assume: M^{z) = —20.29 + ({z — 6). 
We then fit for finding ( = O.SSlo.ii' which results in 
an estimate for M^{z = 10) = —18.0 ± 0.5 mag. 

Alternatively, we assume M^{z) = —19.63 =const. (as 
derived for 2; ~ 8 from our empirical extrapolation), 
a = —1.73, and we fit only for an evolution in the nor- 
malization with redshift relative to the z ^ 8 LF. This 
results in (f)^ = 1.14 x 10"^10"^(^"^^ Mpc"^, with best 
fit T = 0.54l^;?^. Thus, using this extrapolation, the 
normalization of the UV LF from z ~ 10 to z 8 is 
expected to increase by a factor 12. These results are 
summarized in Table [3l 

4.4. Evolution in Luminosity Density From z ^ 10 to 

z^8 

The quantity most easily comparable between observa- 
tions and simulations is the observed luminosity density 
(Pl) above a given limiting magnitude, which is shown 
in Figure [H The LD from the z ^ 10 candidate alone 
amounts to log^Q pl = 24.1^0.7 erg s~^Hz~^Mpc~^. 
However, a more realistic estimate of the luminosity 
density can be obtained from the two possible extrap- 
olations of the UV LF we derived in the previous sec- 
tion, which include the contribution from galaxies at 
^1400 < — 19 mag that are currently undetected. As- 
suming the LF evolution only occurs in the characteris- 
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Fig. 8. — The evolution of the UV luminosity density pL above 
Mi4oo = -18 mag (> 0.06L^=3). The filled circle at z ~ 10.4 is 
the luminosity density directly measured for our only 2; ~ 10 galaxy 
candidate. The two connected dots at 2; = 10 show the range of pos- 
sible LD values, given the two simple, accelerated extrapolations 
of the UV luminosity function described in section 14.31 The red 
line corresponds to the empirical LF evolution from Bouw ens et al.l 
(j2011b|). Its extrapolation to 2; > 8 is shown as dashed r ed line . 
The PL data at z ~ 4 — 8 is taken from Bouwens et al] (|2007l ): 
IBouwens et al. ( 2011b). As can be seen, pL increases by more than 
an order of magnitude in the 170 Myr from z ~ 10 to 2; ~ 8, indi- 
cating that the galaxy population at this luminosity range evolves 
by a factor > 4 more than expected from low redshift extrapola- 
tions. The p redicted pL evolution of the semi-analytical model of 
ILacev et al.l (pOllI^ is s hown as d ashed blue line, and the theoret- 
ical model prediction of ITrenti~et al. (2010) is shown as blue solid 
line. These reproduce the expected luminosity density at 2; ~ 10 
remarkably well. 

tic luminosity, we find for the z ^ 10 luminosity density 
logiQ PL = 24.2 ± 0.5 erg s~^Hz~^Mpc~^, while assum- 
ing the evolution to be driven by a normalization of the 
Schechter function only, we find log^gP^ = ^4.4 ± 0.5 
erg s~^Hz~^Mpc~^. These different estimates are also 
summarized in Table [3l 

Given that the observed luminosity density at z ~ 
8 is \o^,nPL(z 8 ) = 25.6 erg s~^Hz-^Mpc~^ 
(|Bouwens et al.|[2Qllbf ). the inferred increase in luminos- 
ity density in the 170 Myr from z ^ 10 to z ^ 8 amounts 
to more than an order of magnitude. This is a factor > 4 
higher than what would have been inferred from the em- 
pirical relation for the UV LF evolution, which predicts 
an increase by only a factor ~ 3. 

Note however, that such a rapid increase in luminosity 
density is actually predicted by many theoretical mod- 
els. In Figure [8] we also show the luminosity densities at 
^1400 < —18 mag as derive d from the semi-analytical 
model (SAM) of L acev et al.l &)11), and from the theo- 
retical model of iTrenti et al.l (j2010). 

Although there is still some discrepancy on the ex- 
act sh ape of the z > 6 UV LF between the Lacev e t aTl 
(|2011f ) SAM and the observations, the integrated lumi- 
nosity density and its evolution is remarkably well repro- 
duced across the full redshift rang e z ^ 4 — 8 (see also 
discussion in iRaicevic" eTalllMl . The semi-analytic 
model predicts a constant growth in with cosmic time 
at somewhat faster pace than observed, leading to some 
discrepancy between the observations and the model at 
z ^ 7 and z ~ 8, where the observational data shows 
higher luminosity densities. However, the model repro- 
duces our estimates of the 2; ~ 10 LD remarkably well. 

Similar conclusions are reached for the purely the- 
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TABLE 3 

Summary of 2; ~ 10 LF and Luminosity Density Estimates 



logio</>*(^ = 10) 
[Mpc~^mag~^] 


M*(z = 10) 


OL 


[erg s~^Hz~iMpc~^] 


logio P* 
[M0yr-iMpc-3] 


-2.9 (fixed) 
-4.0 ±0.5 

sin| 


-18.0 ±0.5 
-19.63 (fixed) 
^le candidate 


-1.73 (fixed) 
-1.73 (fixed) 


24.2 ±0.5 
24.4 ±0.5 
24 1+°-^ 


-3.7 ±0.5 
-3.5 ±0.5 

q + 0.5 



Integrated down to 0.06L*^3 (M1400 = — 18 mag) 



oretical model of iTrenti et al.l (|2Q1Q[ ). In particular, 
this model predicts the UV LD to evolve at an ac- 
celerated rate at z > 8, being only log;LoPi^ = 24.4 
erg s~^Hz~^Mpc~^ at z ~ 10, in excellent agreement 
with our observed estimates. Since the model is only 
based on the evolution of the underlying dark matter 
mass function, this indicates that an accelerated evolu- 
tion in the galaxy population can be explained even with- 
out the need for a change in the physical mechanisms of 
galaxy formation. 

For furth er theoretical mod el pred i ctions , see also e.g. 
iMao et all (^007); Salvaterra et al.l (120111 ). or iMuhozl 
(|2011[ ). 

It is also interesting to note that the star-formation 
rate densities (see Table [3]) inferred from our data, are 
more than an order of magnitude too low to account 
for the stellar mass densities observed at z ~ 7 in sys- 
tems of similar brightness. With constant star-formation 
over 300 Myr from z ~ lOtoz ~ 7 the observed 
galaxy population would only produce a stellar mass den- 
sity of logio = 4.7 — 5.3 Mq Mpc~^, compared to 
login PA/f(^ = 7) = 6 .6 Mf7^ Mpc~^ as e s timat ed by, e.g., 
iLabbe et al.' ('2010U'£i); Gonzalez et al." (^2010^. If the in- 
ferred mass densities and SFRs are correct, this suggests 
that the majority of the stars found in 2; ~ 7 galaxies 
down to Muy < — 18 mag have to be formed in systems 
below our detection Hmit at z ~ 10 or are younger than 
300 Myr. 

5. SUMMARY AND CONCLUSIONS 

In this paper, we have extended our search for z ~ 10 
galaxies to ~ 160 arcmin^ of public WFC3/IR data ob- 
tained around the GOODS South field (see Figure [1]). 
These data sets have been acquired through the three 
surveys HUDF09, ERS, and CANDELS, and reach to 
varying depths, from Hiqq^ab = 26.9 to Hiqq^ab = 29 A. 
Based on strict optical non-detection requirements and a 
color cut of J125 — i^ieo > 1-2, we search these fields for 
Lyman Break Ji25-dropout galaxies, which are expected 
to lie at z > 9.5. A total of 17 sources satisfy these crite- 
ria. However, 16 out of these sources show strong IRAC 
detections which rule out their being at such very high 
redshifts. Rather, these galaxies are found to have best- 
fit photometric redshifts in the range Zphot = 2 — 4 (see 
section 13.41 and appendix) . They remain undetected in 
the optical due to their evolved stellar populations with 
non-negligible dust obscuration. This shows how impor- 
tant Spitzer IRAC data is for removing contaminating 
lower redshift galaxies. 

Interestingly, these contaminants are essentially only 
detected with magnitudes in the range Hiqq^ab = 24 — 26 
mag. This is ~ 1 mag brighter than the detection limits 
of our bright surveys, which suggests that such evolved 



and dusty galaxies follow a peaked luminosity function 
at these wavelengths. If confirmed by future wide-area 
WFC3/IR data sets, this would indicate that such ex- 
tremely red systems are not as much of a problem for 
z > 9 searches at fainter levels as has been expected to 
date. However, the existence of such galaxy populations 
will make it challenging to use la rge- area WFC3 su rveys 
such as pure- parallel fields (e.g. iTrenti et al.l Mill ) for 
constraining the bright end of the 2; ^ 10 UV LF without 
additional follow-up observations to validate the candi- 
dates. 

Even with our expanded search area, the only > 5cr de- 
tected galaxy with a color limit of i^ieo — [3.6] < 1.6 {2a) 
and thus the only possible z ^ 10 galaxy candidate is the 
same so urce that we reported already in iBouwens et al] 
(|2011af ) . In appendix [B] we additionally note one other 
lower S/N candidate that is suggestive of being at com- 
parable redshift, but requires confirmation from further 
HST WFC3/IR data. 

Interestingly, we would have expected to detect six 
z ^ 10 galaxies in our data, if the UV LF evolved to 
z ~ 10 as expected from lower redshift trends (see §4.2p . 
Thus, the galaxy population appears to evolve at an ac- 
celerated rate beyond z > 8. We infer that the UV 
luminosity density increases by more than an order of 
magnitude in only 170 Myr from 2:~10to2;~8, and 
we are thus likely witnessing the first massive build-up of 
the galaxy population at these early epochs in the reion- 
ization era. The fact that theoretical models based on 
the evolution of dark matter halos do, in fact, predict 
such an accelerated increase in the LD, indicates that 
these rapid changes are mainly driven by an accelerated 
evolution of the underlying dark matter mass function, 
rather than due to a change in star -formation proper- 
ties of these early galaxies (see e.g., ITrenti et al.l 120101 
iLacev et al. 2011). 

The accelerated evolution of the galaxy population also 
has interesting consequences for cosmic reionization by 
galaxies brighter than Muy = —18 mag, the current 
detection limits. With such a sharp decrease in the lu- 
minosity density above z ~ 8, it is impossible for such 
bright galaxies alone to create a reionization history in 
agreement with the high optical depth measurement of 
WMAP and the vast majority of the io nizing fiux has 
to be created by fainter galaxies (see also IBouwens et ahl 

Unfortunately, our knowledge of galaxies at 2; > 8 
still remains very uncertain. However, WFC3/IR offers 
unique opportunities to make significant progress in ex- 
panding the number of galaxies at z ~ 9 — 10 and ad- 
dressing some of the key issues related to early galaxy 
formation and its impact on reionization, even before 
the advent of JWST. In Figure [H we show the depth re- 
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Fig. 9. — Depth required to detect ten z ^ 9 galaxies (blue shaded 
region and dashed line) and ten z ~ 10 galaxies (red shaded region 
and dashed line) as a function of survey area. These regions show 
the expected number of sources calculated using our LF extrapola- 
tions to z ~ 10 from section [4. 31 i.e. assuming only evolution in M* 
or only in to reproduce the detection of our one 2; ~ 10 candi- 
date. The dashed lines show the same, but assuming the standard 
lower-redshift extrapolation of the UV LF (Bouwens et al. 2011b). 
In particular, at 2; ~ 10 the estimated depths differ significantly. 
The range of L*(2; = 10) values of these extrapolations is indicated 
with an errorbar in the lower left. As can be seen, surveys reach- 
ing significantly deeper than ~ 28 mag will be required to detect a 
significant population of z > 8 galaxies. The vertical dashed lines 
indicate the area of one and 50 WFC3/IR fields, respectively. A 
survey with only one pointing would need to reach to ~ 31 mag in 
HiQo^AB to significantly constrain the 2; ~ 10 galaxy population, 
which is out of reach with HST WFC3/IR. Therefore, multiple 
fields are favorable for searching for z ~ 10 galaxies. Additionally, 
comparably deep optical data (> 30 mag) and deep IRAC imaging 
would be required over such fields in order to robustly exclude low- 
redshift contaminants. The identification of 2; ~ 9 galaxies would 
benefit from imaging in different filters (e.g., F140W) than adopted 
in current deep WFC3/IR fields. 



quired to detect ten z ^ 9 and ten z ^ 10 galaxies for 
a given survey area. Such future data sets will have to 
reach significantly fainter than ~ 28 mag to accomplish 
this goal, even for large surveyed areas of 50 WFC3/IR 
fields. This is deeper than currently planned wide-area 
WFC3/IR data (including MCT programs). In fact, to 
characterize the luminosity function at luminosities be- 
low and to set better constraints for reionization, sur- 
veys to fainter than 29 AB mag are really needed. Note 
that with only one WFC3/IR field, a magnitude of ~ 31 
has to be reached to detect a significant 2: ~ 10 popula- 
tion, which is fainter than is practical with HST. 

Based on the current LF constraints, we find that imag- 
ing a few fields provides the best chance to improve on 
current 2; ~ 10 constraints. Furthermore, in order to fur- 
ther constrain the possible accelerated evolution of the 
UV LF, the z ^ 9 regime offers the best opportunity. 
A sizable population of 2; ~ 9 galaxies is expected to be 
seen already down to ~ 29 mag over multiple WFC3/IR 
fields, which can be achieved with the efficient F140W 
filter. Thus, it is likely that already with WFC3/IR, we 
can soon push the frontier of statistical galaxy samples 
with WFC3/IR from z - 8 another - 100 - 170 Myr 
back out to z ~ 9 — 10. 
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APPENDIX 

A. THE INTERMEDIATE REDSHIFT CONTAMINANTS 

Here we briefly summarize the properties of the galaxies that formally met the Ji25-dropout criteria based on the 
HST data alone, but which were found to exhibit very red i^ieo to IRAC colors. This makes it extremely unlikely that 
these galaxies are at z > 9. A summary of these sources is listed in Table IH 

Note that all these galaxies are obviously clearly detected in the IRAC images, with the exception of one source (jD- 
2080646581), which is close to a quintet of IRAC bright galaxies. However, after subtraction of the contaminating flux 
from its neighbors based on a careful convolution of the Hiqq image to the IRAC PSFs (see, e.g., Labbe et al. 2010b!I3; 
[Gonzalez et al.ll2010f ) it also shows a clear detection (see Fig. [TT]) . The best estimate for its color is i^ieo — [3.6] = 
1.6 =b 0.5, which is too red for a likely z > 9 source. 

Interestingly, s i x of these 16 sources a r e pre sent in catalogs of potential passive high redshift galaxies of 
iRodighiero et al.l ()2007[ ) and I Wiklind et al.l ()2008f ). One of these is also detected in X-ray emission and classifled 
as a so-called EXO source hosting a potentially highly obsc ured AGN, particularly if it is an evolved galaxy at inter- 
mediate redshift, as pointed out bv 'Koeke moer et aTl (|2004) . 

We have derived photometric redshifts and mass estimat es for these galaxie s by complementing our HST photometry 
with the IRA C fluxes from the GO ODS MUSIC catalog (jSantini et al.ll2009[ ) or from the SIMPLE images. Based on 
SED flts with iBruzual fc Charlotl (|2003 ) models, these galaxies are conflrmed to be evolved, reddened systems with 
stellar masses around 10^^ Mq at redshifts 2: ~ 2 — 3.5. Note that with the exception of four sources these are all 
detected individually in MIPS 24/im data. 

In Figure [To] we additionally show the surface density of these evolved intermediate redshift sources. Even at their 
peak surface density, they are only found at 0.06 sources per magnitude per arcmin^. Therefore, only ^ 0.4 sources 
are expected per WFC3/IR pointing, providing an explanation for why none of these bright sources are found in the 
three deep HUDF09 flelds. The dearth of such galaxies faintward of i^i60,AB = 26, at ~ 1 mag brighter than the 
completeness limit of our bright surveys, suggests that such extremely red sources become even less frequent (i.e., that 
their LF may have peaked at these magnitudes) and so they could be less problematic for contamination of fainter 




160,AB 



Fig. 10. — The surface density of contaminating sources at redshifts 2; ~ 2 — 4. The surface densities are extremely low. Only ~0.4 sources 
are expected per WFC3/IR pointing. The rapid disappearance of such sources beyond i^i60,AS ^ 26, almost one mag brighter than the 
detection limits of our shallow fields from ERS and CANDELS, suggests that such sources may have a peaked LF at the wavelengths of 
interest in these surveys and so could be less of a problem for contamination for fainter 2; ~ 10 galaxy samples. 
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TABLE 4 

Lower Redshift Contaminants that satisfy the J-Dropout criteria, but show strong IRAC detections. 



ID 



^160 -^125 — ^160 



'^1/21 



MUSIC-IDt references 











ERS 










jD-2162843432 


03:32:16.28 


-27:43:43.2 


23.6 ±0.1 


1.32 ±0.11 


14.6 


0.65 


70081 


_ 


jD-2188742241 


03:32:18.87 


-27:42:24.1 


25.5 ±0.1 


1.36 ±0.23 


12.2 


0.29 


70040 


- 


jD-2226644214 


03:32:22.66 


-27:44:21.4 


25.5 ±0.1 


1.80 ±0.38 


10.9 


0.26 


70104 


— 








CANDELS-Deep 










jD-2487849357 


03:32:48.78 


-27:49:35.7 


25.2 ±0.1 


1.36 ±0.25 


11.4 


0.33 


70314 




jD-2532547516 


03:32:53.25 


-27:47:51.6 


25.3 ±0.2 


1.55 ±0.32 


11.2 


0.26 


70236 


4 


jD-2304648166 


03:32:30.46 


-27:48:16.6 


24.8 ±0.1 


1.23 ±0.11 


29.4 


0.23 


70258 


1 


jD-2387448399 


03:32:38.74 


-27:48:39.9 


24.9 ±0.1 


1.41 ±0.22 


14.9 


0.33 


70273 


1,2,4 


jD-2412344008 


03:32:41.23 


-27:44:00.8 


25.8 ±0.2 


1.43 ±0.30 


11.0 


0.22 


70092 




jD-2158349541 


03:32:15.83 


-27:49:54.1 


24.4 ±0.1 


1.60 ±0.15 


19.1 


0.38 


70316 


1 


jD-2249748085 


03:32:24.97 


-27:48:08.5 


25.7 ±0.2 


1.61 ±0.43 


7.5 


0.33 


70252 




jD-2080646581 


03:32:08.06 


-27:46:58.1 


26.5 ±0.3 


2.01 ±0.67 


6.1 


0.16 












CANDELS-Wide 










jD-2489152264 


03:32:48.91 


-27:52:26.4 


25.2 ±0.1 


1.49 ±0.27 


12.3 


0.26 


70442 




jD-2358952367 


03:32:35.89 


-27:52:36.7 


25.1 ±0.2 


2.27 ±0.94 


7.2 


0.35 


70455 


1 


jD-2351353198 


03:32:35.13 


-27:53:19.8 


25.4 ±0.2 


1.61 ±0.43 


5.8 


0.35 


70484 




jD-2331152057 


03:32:33.11 


-27:52:05.7 


26.2 ±0.3 


1.24 ±0.54 


5.8 


0.20 


70429 


1,3 


jD-2211356269 


03:32:21.13 


-27:56:26.9 


24.4 ±0.1 


1.47 ±0.22 


7.9 


0.44 







Ref erences. — (1) iRodighiero et al.l (I20Q7D : (2) IMobasher et al.l (I2QQ5I) : 
((2Q08h 

t From lSantini et al.l (120091) . 



(3) iKoekemoer et al.l (|20Q1) : (4) IWiklind etlol 



z > 9 dropout searches. This wih have to be confirmed, however, with future wide-area WFC3/IR data, searching 
explicitly for this type of galaxies. 

Note that our finding of a disappearance of passive intermediate redshift galaxies to fainter magnitudes is in agreement 
with Stutz et al.i (;2QQ8i) . who found an absence of passive, red galaxies at z ~ 1.5 — 3 at masses below ~ 10^^ Mq. 

B. POTENTIAL FUTURE EXTENSION OF THE HUDF Z - 10 CANDIDATE SAMPLE 

As pointed out in section O due to the apparent paucity of such sources, it will be difficult to extend the z ~ 10 
candidate lists to large numbers of reliable sources in the near future with HST. In cycle 18, an additional 128 orbit 
imaging survey was granted (PI: Ellis) to further study the z ^1 galaxy population over the HUDF data. However, 
given that only one single field will be targeted and that only a relatively small amount of i^ieo data is currently planned 
to be taken (22 new orbits, relative to the ~58 already available from the HUDF09 and CANDELS programs), we 
can forecast what this program will bring in terms of z ~ 10 galaxy science. Essentially, every 5<j source that will be 
obtained as a result of this new program should already be present as a 4.2cr source in the current i^ieo band data. We 
therefore systematically searched for lower significance sources (down to Act detections) that satisfy our z ~ 10 dropout 
criteria. Essentially, the only additional candidate that we found is already detected at 4.8cr in the current data, and 
is located at RA=03:32:43.01, DEC=-27:46:53.3 (see also Yan et al.l 1201 OD . Note, however, that the i775-band data 
shows very faint positive fiux (~ 1.5cr) exactly at the location of this source, which increases the chance that it is at 
lower redshift. This demonstrates the value of even deeper optical data over these fields for more robust high redshift 
galaxy selections. Given the scarcity of low significance candidates in the HUDF, it is clear, however, that it will be 
very important to image at least one additional field to similar depth in order to further constrain the accelerated 
evolution in the cosmic star- formation rate density (see Fig. [9]). 
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Fig. 11. — Images of all optical undetected sources with J125 — Hiqq > 1.2 mag, but with strong IRAC detections. This rules out the 
possibility that these sources are at 2; > 9. The images are 7.5 arcsec on a side, and show from left to right (1) a stack of B435, Veoe, ^775, 
^850, (2) J125, (3) HiQo, (4) Spitzer IRAC [3.6], and (5) IRAC [8.0]. All these sources are clearly detected in all IRAC bands, including 
[8.0]. Thus, their Hiqq — [8.0] colors are greater than ~ 2, which is very different from the relatively flat color expected for a real 2; ~ 10 
source. Note that IRAC flux of the source jD-2080646581 is heavily contaminated by neighboring sources. In the [3.6] band, we therefore 
show the cleaned image, after subtraction of the contaminating flux, which reveals its clear detection. 



